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According to Anderson’s orthogonality catastrophe, the overlap of the A^-particle 
ground states of a free Fermi gas with and without an (electric) potential decays in 
the thermodynamic limit. For the finite one-dimensional system various boundary 
conditions are employed. Unlike the usual setup the perturbation is introduced by a 
magnetic (vector) potential. Although such a magnetic field can be gauged away in 
one spatial dimension there is a significant and interesting effect on the overlap caused 
by the phases. We study the leading asymptotics of the overlap of the two ground 
states and the two-term asymptotics of the difference of the ground-state energies. In 
the case of periodic boundary conditions our main result on the overlap is based upon 
a well-known asymptotic expansion by Fisher and Hartwig on Toeplitz determinants 
with a discontinuous symbol. In the case of Dirichlet boundary conditions no such 
result is known to us and we only provide an upper bound on the overlap, presumably 
of the right asymptotic order. 

Keywords: Fermi gas, Toeplitz matrix, Hilbert matrix, Szego limit theorem, Fisher- 
Hartwig conjecture 


1 Introduction 

In 1967, P.W. Anderson [l[ discovered that in the limit TV —oo the ground state 
of a free gas of N non-interacting fermions is orthogonal to the ground state of the 
same system but perturbed by an external potential. The asymptotic behaviour of 
the overlap of the two ground states is of the order iV“^ with a constant 7 > 0 
which depends on scattering parameters of the potential. 30 years later, I. Affleck 
Q related Anderson’s orthogonality catastrophe (AOC) to an asymptotic expansion 
of the difference of the ground-state energies of the free and the perturbed system. 
Assuming that the latter is of the asymptotic form cq -l-ci/A^-|-o(l/A^) Affleck argued 
that Cl = 7. Here, the leading term, cq, is called Fumi term and the coefficient, ci, 
of the next-to-leading term is called finite-size energy. 
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In this paper, we study AOC and the difference of the ground-state energies for a 
one-dimensional system which is not perturbed by an external potential but rather by 
an external magnetic field given by a magnetic (vector) potential a. We confine the 
particles to the interval [—L, L], impose periodic or Dirichlet boundary conditions, 
and consider the thermodynamic limit N, L ^ oo where the particle density ^ 
tends to a fixed finite value p > 0. In [^, quantum fluctuations of the current in a 
metallic loop were studied. This led to a logarithmically divergent term similar to 
the orthogonality catastrophe in the present paper. 

It is well-known that in one spatial dimension the Hamiltonian with magnetic field 
can be mapped to a free Hamiltonian (with the same or slightly modified eigenvalues 
as the free problem) by a gauge transformation. Note that while the ground-state 
energy is changed little (or not at all) by the magnetic field, the effect on the overlap 
of the ground states is surprisingly strong. Therefore, we cannot expect that Affleck’s 
identification of the AOC exponent with the finite-size energy holds. 

The main part of this paper is dedicated to the asymptotic analysis of the overlap 
of the ground state of the free Hamiltonian and the one with a magnetic field. For 
periodic boundary conditions, the leading order of the asymptotic expansion of the 
overlap of the two ground states is cN~‘^^ with a finite constant c > 0. A precise 
formulation of this result is given in Theorem 12.31 Here, 6 is half of the integral 
of a (except for additive multiples of tt). This behaviour is remindful of the one 
for potentials. In Theorem 12.51 we give an asymptotic upper bound of this overlap, 
namely with a finite constant c > 0. This second result holds for 

periodic as well as Dirichlet boundary conditions. It is noteworthy that the obtained 
exponent does not depend on the specific form of the magnetic potential a but only 
on its integral 6 which is related to the magnetic flux. 

For the proof of both theorems we rewrite the inner product of the two ground states 
as a Toeplitz determinant with an L-dependent symbol. Due to this L-dependence, 
classical Szegd limit theorems including those for Fisher-Hartwig symbols cannot be 
applied directly. The crucial step of our proofs is that this determinant is factored 
into a determinant, which is asymptotically of order I, times a Toeplitz determinant 
with a discontinuous L-dependent symbol. This discontinuous symbol in turn can 
be interpreted as corresponding to an effective flux which looks like the one induced 
by Dirac’s d-function as a magnetic potential. In the case of periodic boundary con¬ 
ditions the leading 7V-asymptotics of this particular Toeplitz determinant with dis¬ 
continuous symbol is given by a Toeplitz determinant with an L-independent symbol 
of the Fisher-Hartwig class. This exact Toeplitz determinant has been analysed by 
Fisher and Hartwig in their 1968 paper, which they used to bring forward a more 
general conjecture on the asymptotics for a certain kind of discontinuous and singu¬ 
lar symbols Q. Using their result we find the leading asymptotics of the overlap of 
the respective ground states. An upper bound is determined by using the inequality 
det(l -I- A) < where the exponent on the right-hand side is the so-called Ander¬ 
son integral. We study this Anderson integral to obtain the asymptotic upper bound 
for periodic boundary conditions. For Dirichlet boundary conditions we choose a dif¬ 
ferent approach which requires the asymptotic analysis of a determinant of a Toeplitz 
plus a Hilbert matrix. 

In the following section we specify the quantum mechanical model and fix the no¬ 
tation. As a first simple result we obtain the Fumi term and the finite-size energy in 
Subsection o The precise formulation of our main results on Anderson’s orthogo¬ 
nality catastrophe. Theorems 12.31 and 12.51 is given in Subsection 12.31 In Section [3] we 
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prove these two theorems. For this purpose we introduce the notion of generalised 
Toeplitz matrices. In Lemma ld.ll we give a decomposition of the overlap of the ground 
states into a bounded part and a Toeplitz determinant with discontinuous symbol as 
mentioned above. This lemma is proved in Section |4l In Subsection 13.11 we show 
Theorem 12.31 Subsections 13.21 is devoted to the proof of Theorem 12.51 

Finally, we comment on the mathematical work in this field all of which is rather 
recent. The papers focus on AOC (caused by potentials and not magnetic 

fields) and prove imper bounds on the overlap which are all in agreement with Ander¬ 
son’s bounds. In Q and 0, the finite-size energy and its relation to the coefficient 7 
in AOC is studied in a one-dimensional Fermi system. Yet, we do not know whether 
Anderson’s prediction about the decay rate 7 of the overlap is sharp. However, in 
the model of this paper we are able to prove the sharp decay rate. But we also notice 
that Affleck’s correspondence between the coefficient 7 and the finite-size correction 
does not hold here. This should not diminish the value of this fascinating relation in 
the case of a potential. 

2 Model and results 

2.1 The Hamiltonian for a single particle 

We consider a single, spinless particle confined to the interval [—A, L\ C R. The 
Hilbert space is T-Ll ^^([~A,L]), equipped with the usual inner product given 
by {/,(?) := f{x) g{x) dx. We study the system for periodic and for Dirichlet 
boundary conditions where arguments and results for the latter can be taken over 
to Neumann boundary conditions (or more general Robin boundary conditions). We 
denote the eigenfunctions, eigenvalues etc. for different boundary conditions by the 
same respective symbols. It will become clear from the context which boundary 
conditions are used. 

Periodic boundary conditions The Hamiltonian of a single free particle on this 

,2 , . . 

interval is ■= —372 with periodic boundary conditions. The energy eigenvalues 
denoted by Al and the corresponding eigenstates are determined by 

- = ^Wbix) , X G ]-A, L[ , (1) 

<PL{-L)-^dL)=d, (2) 

Then we have, normalised to one in £^-norm, ipLj{x) = ( 2 A )“2 exp (—i^x) and 

Al,j = (^) , j G Z. Note that all eigenvalues except for = 0 are twofold 
degenerate and the eigenvalue 0 is non-degenerate. 

Let a G IT^’^(IR), the Sobolev space of integrable and continuously differentiable 
functions with integrable derivatives. We denote the multiplication operator given 
by (a/)(x) = a{x)f{x) by a as well. In order to have a well-defined operator on 
T-Ll satisfying periodic boundary conditions, we additionally assume a{L) = a{—L). 
The Hamiltonian of the system with the external magnetic (vector) potential a is 
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Ha.L ■= (~i^ —a)^. The spectrum <7{Ha,L) is determined by the eigenvalue problem 



{x) ] iPl{x) = ULtpL^x) , X e ]-L, L[ 


iphi-L) - tpL{L) = 0 , 


d^L 

dx 


i-L) 


d-ipL 

dx 


(L)=0. 


( 3 ) 

( 4 ) 


This eigenvalue problem can be solved using the gauge transformation e'^^ where 
is given by 


$L(a;) 



a{y) dy 


i{y)dy, X e [-L,L] 


( 5 ) 


The quantity $l(x) is called magnetic flux and represents the phase shift due to the 
magnetic field. Since $^(x) = a(x) and a(L) = a(—L), the corresponding eigenval¬ 
ues and eigenfunctions are easily determined: y,Lj = tp]^ j[x) = 

{2L)~i exp^i$i(x) — xj, j S Z. On the one hand, all eigenvalues except 

for 0, which is non-degenerate, are degenerate with multiplicity two if is an in¬ 
tegral multiple of tt and, if $l(T) is a half-integral multiple of tt, the eigenvalue 0 
is also degenerate with multiplicity two. On the other hand, in any other case all 
eigenvalues are non-degenerate. 


Dirichlet boundary conditions For Dirichlet boundary conditions, Eqs. (ED and m 
have to be replaced by 


(Pl{-L) = 0 = (Pl{L) and = 0 =, 

respectively. Moreover, we only assume a € not that a{L) = a{—L). 

The free and the magnetically perturbed Hamiltonian have the same spectrum, Al j = 
= (1^) with j G IN. This is easily shown using the gauge transformation e'^'" 
where is defined as above. The corresponding normalised eigenfunctions are given 

by 


■= 


J Vl 


sindJ-x) for even j 
cos(l^x) for odd j 


and ipL,j ■= 


In contrast with the periodic case, all eigenvalues are non-degenerate and 0 is not an 
eigenvalue. 


Remark 2.1. One could regard b ■.= a' as the magnetic field (whatever the physical 
interpretation of a magnetic field in one dimension is). Due to the condition that a is 
integrable, b has to decay to zero at infinity which, for instance, rules out the case of 
a constant magnetic field. Moreover, the integral of b has to vanish. In the periodic 
case the additional restrictions b{—L) = —b{L) and b{x) dx = 0 for any L occur. 

We want to remind the reader that the gauge transformation maps the Hamiltonian 
with magnetic field, Eq. ©, to a free Hamiltonian (with the same or slightly modified 
eigenvalues as the free problem CD)- However, the eigenfunctions depend on this 
gauge. Therefore the effect of the magnetic potential on the overlap of the ground 
states can be expected to be more drastic than the one on the difference of the 
ground-state energies. 
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2.2 N non-interacting fermions and their ground-state energy 

Periodic boundary conditions Now we consider a system of N non-interacting spin¬ 
less fermions on the interval [—L,L] C IR with periodic boundary conditions. The 
corresponding A^-particle Fock space is the A^-fold antisymmetric tensor product of 
the one-particle Hilbert space, ■F^'^ := A^Hl- The energy of this system is repre¬ 
sented by the Hamiltonian Ha,N,L ■= Ha^L- Note that for a = 0 we recover the 
free Hamiltonian. The ground-state energy Ea,N,L as well as corresponding eigen¬ 
functions of these two N fermion system can be calculated from the eigenvalues and 
-functions of the one-particle system as derived in Subsection 12.II The magnetic flux 
at position L is decomposed as 


$l(T) = riLTT + Sl 


( 6 ) 


with riL e Z and G ] — f, f ] ■ With this notation the ground-state energy of Ha,N,L 

is Ea,N,L = ^ .g hl} C Z for 

odd N and {—m — ul, ..., m — ul — 1} C Z for even N. Here and henceforth we set 
m := ^ for even N and m := for odd N. Note that the ground-state energy 
is twofold degenerate in the cases Sl = N € (2IN — 1) and Sl = 0, iV G (2IN) and 
non-degenerate otherwise. Then the difference Ea,N,L — Eo,n,l of the ground-state 

energies is equal to for odd N and to for even N. 

Therefore the sequence {Ea^N.L — Eq^il^l}l limit points and the asymp¬ 

totic behaviour differs for the subsequences with even and odd values of N, respec¬ 
tively. We find 


Ea,N,L — Eq^n^l 


^+o(i) for odd iV 


as N,L ^ oo with ^ p for a fixed limit density p > 0. Here the parameter 
(5 G ] —-I, is given by mr + S = \ a{y) dy where n G Z is chosen appropriately. 
Thus S = limL_>oo Sl- 

Note that the leading term of the asymptotic expansion, which is of order 1 and 
called Fumi term, is zero. However, the next term in this expansion which is of order 
1/iV, the finite-size energy, is, in general, non-zero for finite values of N. 

Remark 2.2. We want to emphasize that the asymptotic terms may depend, in gen¬ 
eral, how we perform the thermodynamic limit, cf. [3;0| for the case with an (electric) 
potential. 

Moreover, we hnd that the ground state of the Hamiltonians Ha^N,L is the Slater 
determinant Tat.l := tpL,k-nL{x) = <pl,/c(x). 


Dirichlet boundary conditions Analogously to the periodic case we define the N- 
particle Hamiltonians E[q^li,l and Ha^N,L with Dirichlet boundary conditions. The 

Af . 2 

corresponding ground-state energies are Eq^li,l = Ea,N,L = X) (ix) ■ Tf^oy are 

.1=1 

independent of a and therefore the difference of the ground-state energies vanishes 
for all particle numbers N. The ground state of the perturbed system is the Slater 
determinant = A^i = A^i ■ 
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2.3 Anderson's orthogonality catastrophe 

We mainly study the asymptotic behaviour of the overlap of the ground states of the 
free and the perturbed system. In the case of periodic boundary conditions we can 
determine the asymptotics of the overlap rather explicitly: 

Theorem 2.3 (Anderson’s orthogonality catastrophe, periodic boundary conditions). 
Let p > 0 and iV G IN. Assume \ya{y)\dy < oo, max{|nL|} < oo, |(5i| < at 
least for a subsequence of L’s, with |(5| < ^. Then there exists a constant c > 0 such 
that the overlap of the two N-fermion ground states satisfies 




2 


25 ^ 

A^"^(c + o(l)) 


( 7 ) 


in the thermodynamic limit N,L ^ oo with ^ ^ P- In particular, the overlap 
vanishes in the considered limit for 5^0. For (5 = 0 the overlap is equal to 1. 

The proof of this theorem is the main task of this paper and is carried out in 
Section |31 For our proof it is crucial that the matrix on the left-hand side of Eq. ([71) 
is a Toeplitz matrix. The leading term of the asymptotics results from a Szego-type 
limit theorem for a Toeplitz determinant with a discontinuous symbol as given by 
Fisher and Hartwig [^. 

Remark 2.4. Levitov, Lee and Lesovik derived a similar result to Eq. © in a time- 
dependent setup [sl, [lOj. They consider a loop, i. e. an interval of length L with peri¬ 
odic boundary conditions, in a fixed external potential. By a short pulse a magnetic 
flux is induced which leads to a time-dependent vector potential in the Hamiltonian. 
The overlap of the ground states of the system long before the pulse and long after¬ 
wards are compared where the thermodynamic limit L —>■ cx) is taken at fixed Fermi 
energy. 

For reasons that will become clear soon, we, moreover, show a certain upper bound 
to the overlap of the ground states for both boundary conditions. The bound is 
obtained with the inequality In(detA) < — tr(l — A) and calculating the leading 
term of the so-called Anderson integral. For Dirichlet boundary conditions the direct 
calculation of the Anderson integral is omitted in favour of another approach using a 
property of a Hilbert matrix appearing in the asymptotic analysis of the determinant. 

Theorem 2.5 (Upper bound for Anderson’s orthogonality catastrophe, periodic and 
Dirichlet boundary conditions). Let p > 0 and G IN. With the same assumptions 
as in Theorem 1^.31 an upper bound for the asymptotics of the overlap is given by 
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< o(l)) 


( 8 ) 


with a constant c > 0 as N,L ^ oo with ^ t p. 

The proof of this theorem is given in Subsection 13.21 
Remark 2.6. Up to the constant the negative exponent 5'^ of the exact asymptotics 

in Eq. (I7|) and the negative exponent sin^((5) of the upper bound in Ineq. ([5]) satisfy 
“Anderson’s rule”. So they are essentially the sine respectively the arc sine of each 
other. Note that the negative exponent 5'^ is sharp and cannot be improved. 
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Remark 2.7. When considering the overlap for Dirichlet boundary conditions, one is, 
in leading order, led to the asymptotic analysis of 


det 




as TV, L —>• oo . 


Here := 


( I is a Hilbert matrix with 77 G IR, —n 4 IN, and Pat 


projection on an fV-dimensional subspace. In order to determine this asymptotics, 
we would need a Szegd limit theorem for the determinant of the sum of the Toeplitz 
matrix 1 and the square of the Hilbert matrix 77^. We are only aware of asymptotic 
results for the sum of a Toeplitz and a Hankel matrix requiring a certain relation 
between the two symbols which, however, do not apply here. It seems reasonable 
that the overlap is independent of boundary conditions and in particular the same as 
for periodic boundary conditions. 

Remark 2.8. If we would allow a magnetic field b with non-vanishing (finite) integral 
then from the start we can only allow Dirichlet boundary conditions. The logarithm 
of the overlap is now to leading order of the form —N |/5|/(27r). It seems that the 
next-to-leading term is of the order In(A^) and that it is of the same form (possibly 
up to a constant) as above if we replace b hy b — f b in the above asymptotics for 
vanishing total magnetic field b. 


3 Asymptotic analysis of Toeplitz determinants 


Let / : [—L, L] —>■ C be a complex-valued (essentially bounded) function. The oper¬ 
ator of multiplication by / is denoted by / as well. A generalised Toeplitz matrix 
with symbol / is defined as the N x Wmatrix 

Tjv(/) := 

where is an orthonormal basis of the underlying Hilbert space. The correspond¬ 

ing determinant det(TAr(/)) is referred to as the (generalised) Toeplitz determinant. 
Here the generalisation is with respect to the choice of the basis. In our context we 
choose the eigenbasis of the corresponding free Hamiltonian. In the case of periodic 
boundary conditions this basis consists of plane waves and the matrix elements are 
the Fourier coefficients of /. In particular, they only depend on the difference j — k of 
the indices and we recover a Toeplitz matrix in the usual sense. For Dirichlet bound¬ 
ary conditions, the emerging matrix is a sum of a usual Toeplitz and a Hankel matrix. 
For a deeper study and an overview of results for Toeplitz and related matrices see, 


for instance, 11 M- 


Hereinafter, we study symbols of the form / = e'®^. In the case of periodic bound¬ 
ary conditions is given by 


gL{x) := ^l{x) - —X 


for X G [—L, L] 


(9) 


with and <5 l as defined in Eqs. and respectively. The function has a 
jump of height 2nLTr at the boundary x = ±L. This becomes particularly simple if we 
take, at least informally, the discontinuous magnetic potential a{x) := 2 ^l{L)S{x). 
More precisely, we set 


grix) := $l(L) sign(a;)- -x for x G [-L,L] 


( 10 ) 
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which is reminiscent of an integral of a. As gi,, gL has a jump discontinuity of height 
2 nL7r at the boundary x = ±L, but also at a; = 0 with height 2 $l(L). 


For Dirichlet boundary conditions we do not have the term —^x. Hence we set 
5 z,(x) := $l(x) and 5 l(x) := $l(L) sign(a:) for x€[-L,L]. (11) 


Here the jump discontinuities at 0 and at the boundary are of the same height 2$^ (L). 


Now we have 


det{{(pL,j,tpL,k)),k=i,...,N = |det(r/v(e'9^))|^ . Note that the sym¬ 


bol of this Toeplitz determinant depends on the parameter L, which also tends to 
infinity, unlike in classical Szegd limit theorems. To circumvent this problem Dn,l '■= 
det (T^v (e‘®^)) is factorised into the Toeplitz determinant ■= det (TAr(e‘®^)) 

with the symbol e'®^ and a bounded term as the main step in the proof of the theo¬ 
rems. This Toeplitz determinant Dn,l in turn can be rewritten as a usual Toeplitz 
determinant with a Fisher-Hartwig symbol times a remainder term of order 1. More 
precisely, we use the following lemma. 


Lemma 3.1. For the system with periodic or Dirichlet boundary conditions assume 
f-oo \y^(.y)\‘^y < max{|ni|} < oo, and |(5i| < c for some 0 < c < ^. Then we 
have 

|det(T^(e‘®^))|' = Cn,l ■ |det(rw(e'®^)) 


for any N and L. In general, Cn,l depends on N and L. Moreover, there are positive 
constants po, (5i, and 62 such that 

N 

0 < (5i < Cn,l < 62 < 00 for all N and L with — < po • 

ZIj 

The proof of this lemma is postponed to Section 0] It is based on a factorisation 
of the determinant and a decomposition of some of the resulting matrices. 

Remark 3.2. For simplicity we henceforth assume that the quotient ^ is kept at 
a fixed finite value p > 0 when we take the thermodynamic limit. However, this 
condition may be relaxed to lim 7 v^L_>oo ^ = P nnd the following statements still hold 
true. 


3.1 Proof of Theorem 12.31 

In the following we show Theorem l2. 31 using the decomposition of the determinant as 
given in Lemma 13.11 

Proof of Theorem \2.tA According to Lemma l3.li the leading asymptotics results from 
the Toeplitz determinant with the discontinuous symbol e‘®^ with g^ given by Eq. m- 
Thus the leading term arises from 


Dn,l 


sin((5L) 


det 




( 12 ) 
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where m = • This can be verified by a straightforward computation using Eq. (|5|) . 

With the abbreviation Sm ■= we obtain the N x iV-matrix 

OL—TTm 



/ So 

S-1 

S-2 

S-N+2 

S-N+l\ 

Tjv(e‘»") = 

Sl 

So 

S-1 

S-N+3 

S-N+2 


\SAr_i 

SN-2 

SN-3 • 

Sl 

So / 


This is a Toeplitz matrix with the symbol e'®^ which is discontinuous at zero and at 
the boundary, but has no singularities. 

If a is compactly supported, Sl is constant for sufficiently large values of L and 
the constant value is S. If we relax the condition that a has compact support, the 
L-dependence of the symbol has to be taken into account. However, with ||Xa||i < oo 
and II a 111 < oo, a has to decay faster than ^ at infinity. Therefore \Sl — S\ = o{y) 
as L —^ oo. Thus the first error term decays faster than L~^. Consequently, we 
may assume the magnetic potential to be compactly supported from now on. In this 
special case there is an L < oo such that Sl is constant (i.e. does not depend on L) 
and is equal to 5 for every L > L. 

By the above arguments, the parameter 5l can be considered to have the constant 
value 5 = limL_>oo for sufficiently large L. The Toeplitz determinant we obtain this 
way has already been considered by Fisher and Hartwig in 1968 Q as an example 
where the conjecture they brought forward in the same paper, now called Fisher- 
Hartwig conjecture, holds. They explicitly calculate the leading asymptotic term in 
the limit N ^ oo using Cauchy’s determinant formula. So we find 

|det(rAr(e’9'^)) 1^ - + o(l)) 


as iV, L —>■ oo with ^ —5> p. Note that the constant E depends on S and is given 
explicitly in Q. Related results on the Fisher-Hartwig conjecture can be found in 
3 M 1 3 ■ Merging the above results and Lemma 13.11 we obtain the assertion of 
Theorem 12.31 □ 


Remark 3.3. In his study on Anderson’s orthogonality catastrophe [l^, P. Lloyd 
constructs a model unitary transformation between the eigenfunctions of the free 
Hamiltonian Hq and the perturbed Hamiltonian Hq + V. The matrix elements E 
Eq. (A.2.4)]of this unitary transformation with respect to the free eigenfunctions are 
of the same form as in (TT^ with 6 ^ replaced by the phase factor of the potential, 
which allowed him to justify the exponent in Anderson’s orthogonality catastrophe. 


3.2 Proof of Theorem 12.51 

Periodic boundary conditions and the Anderson integral We first show Theo¬ 
rem 12.51 for periodic boundary conditions. To this end we study the asymptotic 
behaviour of the corresponding Anderson^integral. Again, it suffices to study the 
asymptotics of the Toeplitz determinant, with the discontinuous symbol e'®^. 

For simplicity we assume ul = 0 and N to be odd in the following and set m := . 

For riL ^ 0 the area of summation has to be shifted by to the left and for even 
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N = 2m the term with j = m has to be removed from the sum. With the two pro- 
jections Pat := X) iV’L.j, ■)‘PL,j and Hat— X ^Plj, ‘Plj have 

j— — m j——m 


Dn,l 


det({(^Lj,e*®"'<^L.fc))^-= det (PArnArPAf) 


The inequality det(A) < exp(—tr(l — A)), which holds for trace class operators A, 
yields 


det^^PAfllArPAr) < exptr (^Pat(1 - nAr)PAr) ) = exp(-lAf,L) 


where Xat^l := X 

j^-m 


oo _ cy —m —1 

k— — oo 

sin(SL) 


is the An¬ 


derson integral. Since {^L,j,e^^^(pL,k) = ; the individual summands of the 

Anderson integral only depend on the difference of the indices j — k. Rescaling and 
renaming the summation variables we obtain 


N oo 


sin {Sl) 


1n,l = 


1 


N oo 


+ EE- 


1 


=1 fe=i {k + j — \ (—k — j + 1 — 

k ,, A 1 


= E 


+ N y — 

k^+i (fc -1 - %) 


N 


+ E 




-N y 


1 


fc=Ar+i 




2 ■ 


On the one hand, the second and the fourth sum on the right-hand side are mono¬ 
tonously decaying at least as ^ for large N. On the other hand, the first and the 
third sum are diverging for N ^ oo and their respective asymptotics are In(A^) + 
0(1). Hence the leading term of the asymptotics of the Anderson integral is In,l = 


_2 „;„2 
TT^ 

the above derivation yields 


^ sin^(i5) In(Af) -|- 0(1) as N, L ^ oo with ^ p > 0. Together with Lemma [3A] 


< o(l)), 


det((XL.j-,V'L.fc)),-fe_„,...,„ 
thus finishing the proof of Theorem 12.51 for periodic boundary conditions. 


Dirichlet boundary conditions and the Hilbert matrix For the Dirichlet case we 
use a different way to prove the respective assertion of Theorem 12.51 This approach 
has the advantage that we have equalities until the very last step where the only 
estimate of this proof comes into play. It shows a possible way to tackle the problem 
of the exact asymptotics of the overlap as already mentioned in Remark 12.71 

For simplicity we consider the particle number to be even and set N =: 2M 
throughout this subsection. For odd particle numbers N some of the following for¬ 
mulae are slightly more complicated, but the general method works as well. In order 
to further simplify ■= detTA^(e‘®^) we have a closer look at the matrix elements 
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e^^^ipk) that appear in this determinant. By the trigonometric addition theorems 
we find 


ipL,j{x)(pL,k{x) = < 


COS 



2t(cOs(^^x) -COsf^ig^X 


for j, k even, 
for j, k odd , 

fT j for j even and k odd , 

^ sin^ ^^^^*^^ x'j + sin^ ^^^^^^ ^ for j odd and k even. 


L _ L _ 

We obtain J cos(^x) dx = 0 for m ^ 0, f cos(^x) dx = 

-L -L 

2Lcos($l(L)) for TO = 0 and sin(^x) dx = i^sin($i(L)) for odd 

TO. Merging these results yields 

"cos($l(L)) forj = fc, 

(iOr e‘®-C.r fo = <! - ^ ^ 

^ sin($L(L)) [j^ - j^] for j odd and k even, and 

^ 0 for J ± A: even. 

In symbolic notation the resulting determinant has the form 


D 


N,L 


det 


det 


/Cl Cl 

C1S2 

C1C3 

C1S2M \ 


2 


S2C1 

S2S2 

S2C3 

S2S2M 




\S2M 

S2MS2 

S2MC3 

S2MS2m) 




/ Cl Cl 

C1C3 


C1C2M-I 

Cl 52 

C1S2M \ 

C2M-1C1 


C2M-IC2M-I C2M-IS2 

C2M-1S2M 

S2C1 



S2C2M-I 

S2S2 

S2S2M 


\ S2 MCi 


S2MC2M-1 


S2MS2 


S2MS2M j 


where we interchanged columns and rows for the last identity. Due to the particular 
structure of the matrix elements we obtain 




, cos($l(L)) 1 -^sin($L(L))iJ];^ 

_|isin(<i>A(L))if'^ cos(<I>i(L))l 


Here the two M x M-matrices and have the entries \ 

and 5 respectively. The determinant is further reduced to 


Dn,l 


det [ cos^(<I> a(L)) 1 H —- sin^($A(T)) 


det( 1-J sin^((5L) Km 
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using the 7 r-periodicity of sin^ and ^l{L) = n^Tr + Sl with nz, G Z, |dz,| < The 


matrix Km is defined by {KM)jk ■= jk X) 


i=M+i 


. Using the ele¬ 


mentary formula 


jk 




1 1 


1 1 

i-^-j *-pi_ 


— k l—-^-\-k 


we 


obtain the decomposition Km = Kj^ + K'^ + Kj^ + K^ . Here the matrices are 

OO OO 

given by (Kj^ ) := 4 Z) i-h-j i-k-k ’ )jk ■“ ~I S i-k+j i-i-k ’ 

l=M+l 2 J 2 J l=M+l = 2 


i^M ) jk '■ 4 Z l-^-j l-^+k {^M ) jk '■ 4 Z i +4 /_ i+fc ■ 

i=M+l 2 J 2^ J l=M+l ^ = 


Intro¬ 


ducing the operators A := 


j,ke]N 


B := 


(^Pi+fe) 


j,fceiN 


and the orthogonal 


M 


projection Pm '■= Z (V’i.fe) ' )VL,k we obtain 

fc=i 


Km~ = iPMA*{t - Pm)APm , K+- = -IPmB{1 - Pm)APm , 

Kit = -\PMA*{t - Pm)BPm , K++ = IPmB{1 - Pm)BPm ■ 

Note that B is self-adjoint, B* = B. In order to estimate the determinant we need 
the trace norms of these operators. To begin with, 


. ^ M OO 

\Kttt = -\\PMB{t-PM)BPM\l = jY. E 


fc=M+i \J ^ ^ 2 ) 

Analogously, we obtain 

- ^ M 00 

II^M-|li = 7 II^Mdl*(l - PmMPmIIi = tE E 


= 0 ( 1 ). 


j — 1 




M 


“Pi: 


1 


,=i k=j (* + 1 ) ^ 


= -ln(iV) + 0(l) 


for large N = 2M. The mixed terms are estimated via the Cauchy-Schwarz inequal¬ 
ity: 


< I \\PMA*{t - PM)h \\PMB{t - PM)h = OiVHm) . 


Since Kj^ gives the leading term for large N, we need the operator norm of K^ . 
To this end we first rewrite the matrix elements as 

1^1 1 1 “ 1 1 


Next we define the operator Qm ■ —>■ by {QMu)j •= UM-j for any j = 

Obviously, &m = 1 In particular, Qm is unitary. 

Furthermore, we need the infinite Hilbert matrix PI 1 := (—I —r 1 • Then 

2 \ 3 +k-^ J j^kdK 


Km = j^mPmH^ iPmQm and IliF 


M 


\ op 


<i 


H 


op 


= ^ by the well-known 
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result 


H 


op 


= TT [13, We conclude that the operator 1 — %siv?{5L)K 


M 


is invertible since '^siT?{5L)K^ ||^p < ^ for 0 < |5 l| < Therefore this operator 
yields the leading term and we have 


det({(pL,j,V’L,fe)),fe=i 


,iV 


det (1- - sm‘^{SL) K 


det 1 — 


t - - sin^(S l)K 


M 


TT^ 

-1 


M 


— sm^{6L) [K++ + K+f + KI+] 


The second determinant on the r. h. s. is bounded (with similar arguments as in 
the proof of the Lemma EH). The asserted upper bound is obtained by using the 
inequality det(l + A) < exp(tr A), 


det( 1- - sin^((5L) K 


M 


< exp 

= exp 


^ sin^{S l)\\Kj^ lli^ 
^^^^ln(iV) + o(^/h^)) , 


which completes the proof of Theorem 12.51 


4 Singling out the effective flux - proof of Lemma 13.11 

In the following we prove the first assertion in Lemma |3.II As a first step, we show 
that TAr(e'®^) has a bounded inverse. This is easily concluded from the following 
general properties of a Toeplitz matrix. Due to the special structure these properties 
of a Toeplitz matrix are taken over from its symbol. We state them without proof. 

Proposition 4.1. (i) The map f T^i^f) is linear. 

(ii) If f is real valued, T^if) is self-adjoint, i. e. T]s[(f)* = Tpf(f). 

(Hi) If f > 0 (f <0), then Tj^^f) is positive (negative) semidefinite. 

(iv) //±Re(/) > 5 > 0, then TN{f) is invertible with ||T/v(/)~^|| < |. 

Remark 4.2. Assertion ifwl) is a corollary of the Lax-Milgram Theorem. 


4.1 Periodic boundary conditions 

The real part of the symbol of Tjv (e'®^) is 


(-1)"^ Re (e'9^ W) = (-1)"^ Re > cos{5l) 

for all X e [—L, L], since |5 l| < ^ and ghix) € ]—ultf — 5l, —riiTr] n [n^Tr, niTr + 5l\- 
Then Proposition 14. llfnil) implies that T 7 v(e'®^) is invertible and 




< 


|cos((5l)| ' 

In order to investigate (e'®^) we rewrite it in the following way: 

rjv(e'®^) = Tiv(e'»^) [l + rw(e'®^)”'(r^(e*9^) - T^(e'®^)) 


( 13 ) 
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Eventually, we want to study determinants and therefore need the trace class prop¬ 
erties of the difference A^r := Tjv(e‘®^) — T 7 v(e‘®'"). Since the trace norm can be 
computed most easily for semidefinite matrices, we first decompose Ajv into hermi- 
tian matrices: A^r = T 7 v(e+) + TM{e~) — iTAr(/+) -|-iTjv(/“)- This decomposition is 
used to obtain an upper bound for the trace norm of A^y: 

IIA^IIi < ||r^(e+)||, + ||T^(e-)||, + ||T^(/+)||i + ||T^(r)||i • 

By linearity, it suffices to decompose the respective symbols. We start by splitting 
the symbol into real and imaginary part 

gi 9 i(a:) _gi 9 i,fo) _ cos(gL( 2 ^)) - cos(ffL(x)) + i [sm(g£(x)) - sm(g£(x))] . (14) 

We recall the trigonometric identities cos(u) — cos(v) = — 2 sin(2i^) sin(^^i^) and 
sin(M)—sin(u) = 2 cos(^^^) sin(^^i^) which are used to rewrite the real and imaginary 
part. Then the difference gL{x) — ghix) and the sum gL{x) -b gL{x) occur. We first 
compute the difference and obtain 

-\-L 

9 l{x) - gL{x) = - j a{y)dy =:-^l{x) for x>0 and (15) 

X 

X 

gL{x)-gL{x) = J a{y)dy =: <^l{x) for x< 0 . (16) 

-L 


For the sum we find 


gL{x)+gL[x) = ^'^{x)—2—x for x>0 and 
gL{.x)+gL{x) =-^J^iyx)-2^x for x<0. 


Rewriting the trigonometric functions in the real part on the r. h. s. of Eq. da we 
obtain 


5l 


cos(( 7 l(x)) — cos( 5 L(a;)) = —2sini -$))(x)-—x ) sin —(x) ) for x>0 and 


1 


L 7 "“7 2 

cos(gL(x)) - cos(gL(a;)) =-2sin^-i$£(x) - ^x^ sin^i$+(a;)J for x< 0 . 

Analogous computations yield for the imaginary part 

sin(gj;,(x)) — sin(gi(x)) = — 2 cos^i$7a:) — ^x^ sin^i$£(x)^ for x>0 and 
sin( 5 L(x)) - sin(gL(x)) = 2 cosQ$))(x)- b ^x j sinfi^^x) j for x<0. 


This yields a decomposition into four hermitian Toeplitz matrices with the symbols 


5l 


e^{x) := 20 (x)sin( -$ 7 x) — j-x ) sin( -<i)^(x) ) , 


1 
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e-{x) := 20 (-a;)sinQ$^(a:) + sinQ$+(a;)^ , 

f+(x) := 20(x)cos^i$J(x) - ^x^ sinQ$^(a;)^ , 
f-{x) := 20(-a;)cos|^i$£(x) + ^x^ sin|^i$+(x)^ , 

{ 0 for X <C 0 j 
1 for a; > 0 . 

However, these four hermitian matrices are indefinite in general and we have to 
estimate their trace norm. In view of Proposition ld.Hfml) we decompose the symbols 
in their positive and their negative part, e. g. 

e>(a;) := max( 0 , e^(a:)) and e^(a;) := min( 0 , e^(x)) . 

Thus we split each of the four hermitian Toeplitz matrices into a positive semi-definite 
and a negative definite Toeplitz matrix. We give the arguments for the decomposition 
of e+j the estimates for the other three cases are obtained similarly. We dehne the 
sets 


/> := {a; G [ 0 ,T] | e'''(a;) > 0 } and /< := {x G [ 0 ,L] | e“'"(a;) < 0 } 

which yields a disjoint partition of the interval [0,L]. Note that e“''(a;) = 0 for all 
a; < 0. Recall the definition of in (HH) and the simple inequalities |sin(a;)| < 1 
and I sin(a;)| < |a;|. Using them we find for the matrix elements 




i> 


L 

J J \a{y)\dydx 

/> X 


sin 



da; 


and analogously 


L, 

\{TNiet))kk\ - ^ J J Hy)\dydx. 


7< X 


Since />U7< = [0,L] and \a{y)\dydx = J^y\a{y)\dy, we have 

L 

,, ,, ,, n N P 

||7w(e+)||^ < TN{e^) ^ + ||Tw(e+)||^ < — J y\a{y)\dy 

0 

Bounds for the matrix elements of TN{e~) and TN(f^) are found likewise: 

L 0 

N f N f 

||7w(/+)||i < ^ J y\a{y)\dy and ||T 7 v(e")||^-7 ||TAr(/-)||^ < — J \ya{y)\dy. 


-L 
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Altogether we have 


L 

IIAatIIi < y / \y^^y)\'^y- 

-L 

Merging the estimates given in Eqs. (US and (El) we obtain the bound 

1 ^ , M , 

-L 


(17) 


In order to ensure that the determinant is non-zero, the right-hand side can be made 
as small as we wish, e. g. if the quotient ^ is considered to be small. We conclude 
that the determinant is estimated by 

0 < <5i < |det(l-tTA,(e'90”^ T < ^2 < oo 

with appropriate constants 5i and 82 which finishes the proof of Lemma l3.II 
Remark 4.3. The bounds on the matrix elements of T]sf{e^) and 7Ar(/^) hold for a 
general orthonormal basis of essentially bounded functions: 


{TN{e^)),k 


< 2||(pj||oo||<Pfe||oo / y |a( 2 /)|dy and 


{TnU < 2||(/7j||oo||<Pfe||oo / |y||o(y)|dy. 


-L 


4.2 Dirichlet boundary conditions 

The proof for Dirichlet boundary conditions is similar to the one shown above, except 
that the term is not present in the Dirichlet case. The difference is that now the 
decomposition 

TNie's^) - rjv(e‘®^) =En- iF+ + iF" 

is into three positive semidefinite Toeplitz matrices, instead of four, with symbols 
e{x) := 20(-a:)sin(i$£(a:)) sin(i$+(x)) , 
f+{x) := 20(a:)cos(i$)-(x)) sin(i4>^(a;)) , 
f-{x) := 20(-a:)cos(i$^(a:))sin(i4>+(a;)) . 

The assertion of the Lemma o then follows by adopting the arguments of the pre¬ 
vious proof. 
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